Reconciling the NRQCD prediction and the J/ift — >■ 37 data 
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Abstract 

It has been a long-standing problem that the rare electromagnetic decay process J ftp — > 37 is 
plagued with both large and negative radiative and relativistic corrections. To date it remains futile 
to make a definite prediction to confront with the branching fraction oiJ/ip — > 37 recently measured 
by the CLEO-c and BESIII Collaborations. In this work, we investigate the joint perturbative and 
relativistic correction (i.e. the 0(a s v 2 ) correction, where v denotes the characteristic velocity of 
the charm quark inside the J/ip) for this decay process, which turns out to be very significant. 
Q-| After incorporating the contribution from this new ingredient, with the reasonable choice of the 

input parameters, we are able to account for the measured decay rates in a satisfactory degree. 
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Calculating the ortho-positronium (o-Ps) annihilation decay into three photons, together 
with the hyperfine splitting of the positronium, has played a pivotal role in developing the 
bound-state QED formalism, particularly in establishing the modern theoretical framework 
such as nonrelativisitc QED (NRQED) [lj. After incorporating a number of higher order 
corrections [2|-0] , the theoretical prediction to the decay rate for o-Ps— » 37 is consistent with 
the measurement [EJ to an impressive accuracy, which can be viewed as a great triumph of 
NRQED. 

The QCD analogue of the o-Ps, the J ftp meson, has also occupied a special stage in 
the making of the Standard Model. The numerous decay channels of the J/ip have been 
intensively studied in the past four decades, among which the rare decay J/ip — > 37 is 
quite interesting. Despite considerable efforts, experimentalists have not observed this decay 
channel until quite recently. In 2008, the CLEO-c collaboration measured the branching 
fraction of this process for the first time 0] and obtained Br(J/-0 — > 37) = (1.2 ± 0.3 ± 
0.2) x 10 -5 . Very recently, with much large J/ip samples, BESIII Collaboration refined the 
earlier measurement and gave Br(J/V> 37) = (11.3±1.8±2.0)xl0- 6 [7]. Both experiments 
are compatible with each other, nailing down the branching fraction for this rare decay to 
be about 10~ 5 . 

The theoretical investigation on J/ip — > 37 preceded the experimental discovery long 
earlier. By the early 1980s, both the first-order radiative and relative corrections have al- 
ready been calculated in the context of the potential model [1 T3] . Just analogous to the 



very successful application of NRQED to positronium, nowadays the nonrelativistic QCD 
(NRQCD) effective field theory [lj becomes the standard tool to analyze the quarkonium 
spectrum, decay and production [ll|. In particular, the influential NRQCD factorization 
approach [l2| has superseded the phenomenological potential model as a modern and sys- 
tematic framework. To the best of our knowledge, up to the error of relative order v 4 , the 
decay rate of the process J/ip — > 37 in the NRQCD formalism can be written as 

r(j/^3 7 ) = 8(?r2 ~ 9) / ca3 1 (olx f g • JM I 2 { 1 ~ 12 -630 - 
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where e c = | and m c represent the electric charge and the mass of the charm quark. (0|x^cr • 
e*ip\J/ifj(e)) is the lowest-order (LO) NRQCD J/ip-to- vacuum matrix element, where e 
represents the polarization vector of J/ip. Up to an error of relative order v 2 , this matrix 

element can be approximated by y^i?j/^,(0), where Rj/^(0) denotes the wave function at 

the origin for the J/ip in the potential model. A useful ingredient of this NRQCD formula is 
the inclusion of the leading relativistic correction, characterized by the ratio of the following 
NRQCD matrix elements: 

= (QixW(-ffi)ViJM m 

[ ' J/lp ~ m 2 c (0\x^cr-e*^\J/^) ' K ' 

where T? stands for the left-right asymmetric covariant derivative. The Gremm-Kapustin 
relation [l3j, which stems from the NRQCD equation of motion, can be employed to express 
(v 2 )j/^ in terms of the J/tp mass and the charm quark pole mass. Empirically, one may 
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FIG. 1: The renormalization scale dependence of the branching fraction for J/tp — > 3j. The 
LO prediction is represented by the black solid line, the LO result plus the 0{v 2 ) correction 
by the dashed line, the LO result plus the 0{a s ) correction by the dotted curve, the prediction 
incorporating both the 0(a s ) and 0(v 2 ) corrections by the dot-dashed curve. Finally, the complete 
NRQCD prediction incorporating the new 0(a s v 2 ) correction is represented by the red solid curve. 
For the sake of comparison, the BESIII measurement of the decay branching fraction is represented 
by the green band, with statistical and systematic errors added in quadrature. 



expect that {v 2 )j/^ ~ 0.3, compatible with the typical velocity of charmonium deduced from 
the potential model. 

The leading term in (JT|) is translated from the analogous formula for the ortho-positronium 
decay to 3 photons, originally derived by Ore and Powell in 1949 [lij • The 0(a s ) correc- 
tion [10] was adapted from the 0(a) correction to o-Ps— > 37, first correctly calculated by 
Caswell, Lepage and Sapirstein in 1977 [2], or from the 0(a s ) correction to J/if> — > 3g, 
first calculated by Mackenzie and Lepage in 1981 [sj. Note at that time the 0(a s ) short- 
distance coefficient can be obtained only with very limited precision. The relative order- 
v 2 correction was originally obtained by Keung and Muzinich in 1982 |9]. With the aid 
of the Gremm-Kapustin relation, their result is consistent with the value given in ([I]), 

iPf^ 2 }./^ -5.32<^ 0. 

Our current understanding of J/tp — > 37 is essentially not much better than three decades 
ago, since to date the 0(a s ) and 0(v 2 ) corrections remain to be the only known corrections. 
Upon confronting ([I]) with the J/ip — > 37 data, the agreement seems far less satisfactory 
than that for its ortho-positronium cousin. The situation is most clearly illustrated in 
Fig. [TJ With some reasonable choices of the NRQCD matrix elements, the LO prediction 
is several times greater than the new BESIII measurement, though capturing the correct 
order of the magnitude. After including the 0(a s ) correction, the decay rate becomes 
significantly lower than the measured one, and even turns negative for most plausible range 
of the renormalization scale entering the a s . The situation deteriorates desperately after the 
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FIG. 2: Some representative C?(a s ) Feynman diagrams for J/ip —> 3j. They include the a) tree 
diagram, 6) self-energy, c) outer vertex, d) inner vertex, e) double vertex, /) ladder, and </) wave 
function renormalization. 



(9(w 2 ) correction is further included, in which for virtually all the plausible input parameters, 
the decay rate becomes deeply negative thus loses physical significance. 

The utter failure of predicting the decay rate for J /if) — > 37 may imply the breakdown of 
the NRQCD for this decay process. The symptom seems to be rooted in the uncomfortably 
large negative radiative and relativistic corrections. In 1997, Braaten and Chen estimated 
the asymptotic behavior of the higher order radiative corrections to the coefficient associated 
with the LO NRQCD matrix element, and found the series are badly divergent due to the 



large residue of the u — — 5 infrared renormalon 16 



Since the large negative relativistic correction appears to be even more troublesome than 
the radiative correction, in this work we attempt to investigate the 0(a s v 2 ) effect, by com- 
puting the first-order radiative correction to the short-distance coefficient associated with 
( v2 )j/i>- Thus far, the 0(a s v 2 ) correction has only been available for a few reactions in- 
volving S'-wave quarkonium, e.g., J/if) — > e + e~ 17, [lij], B c — > W [19J, rj c — > 77 20, 21], 
rj c — > light hadrons j2l|, 22| . e + e~ — > J/ip + r] c |23|. The 0(a s v 2 ) correction turns out to be 
insignificant for most of the aforementioned processes. Nevertheless, for the hadronic decay 
of r] c , the 0(a s v 2 ) correction is found to be important [21]. Therefore, it is also of some 
interest to assess the impact of the 0(a s v 2 ) correction for J /if) — > 37. 

The 0(a s v 2 ) coefficient is represented by the entity in the parenthesis in (JT|). The first 

term that depends on the NRQCD factorization scale /!/, -jPln^-, can be inferred solely 
based on the knowledge of the one-loop UV divergence of the lowest-order NRQCD matrix 



element [T7J and the Ore-Powell result. Therefore, our central task is to calculate the 
unknown constant G. 

We proceed to briefly describe our calculation of the desired NRQCD short-distance 
coefficients. The technical details will be presented in a long write-up. Rather than invoking 
the standard perturbative matching approach, we resort to a shortcut by directly extracting 
the hard region contribution of a loop diagram in the spirit of method of region [2~~ 
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For the process at hand, the NRQCD factorization also holds at the amplitude level. We 
start by considering the on-shell amplitude for cc — » 37 through 0(a s ). Some typical Feyn- 
man diagrams are shown in Fig. [2j At the tree level, there are 6 diagrams for this process. 
At next-to- leading order in a s , there are 12 self-energy diagrams, 12 outer vertex diagrams, 
6 inner vertex diagrams, 12 double vertex diagrams, together with 6 ladder diagrams. The 
Mathematica package FeynArts 25[ is utilized to generate the Feynman diagrams and 
the corresponding amplitude. 

We use Dimensional Regularization to regularize both UV and IR divergences. It is 
convenient to employ the covariant spin projection technique j26[ to enforce the cc pair to 
be in the color-singlet and spin-triplet state. The color and Dirac trace algebra are handled 
by the package FeynCalc j27[. Prior to carrying out the loop integration, we expand 
the amplitude in powers of relative momentum between c and c through the second order, 
followed by the projection of the S'-wave orbital angular momentum state. 



We then employ two different versions of the self-written Mathematica codes [28| to 
carry out the partial fractions, to reduce the general higher-point tensor one-loop integrals 
into a minimal set of scalar integrals up to the four-point. Up to this step, some scalar 
integrals may contain propagator of cubic power. The package FIRE 29J is then employed 
to reduce these unconventional scalar integrals to the standard ones. All the required one- 
loop master integrals can be found in Ref. 0]. 

After summing the contributions from all the diagrams, and taking into account the 
wave function renormalization of the charm quark, the ultimate QCD amplitude becomes 
completely UV finite, nevertheless contains a piece of unremoved logarithmic IR divergence 
at the relative order a s v 2 . The occurrence of the IR divergences in the hard region at this 
order is just as expected, physically because of the breakdown of the color transparency once 
beyond the leading order in v. This IR divergence can be absorbed in the MS-renormalized 
NRQCD matrix element, and consequently, the 0(a s v 2 ) short- distance coefficient now ex- 
plicitly depends on the NRQCD factorization scale /i/, with a natural range from m c v to 
m c . 

Since the differential short- distance coefficients inferred at the amplitude level are IR 
finite, one can safely return to 4 spacetime dimensions. Squaring the amplitude, summing 
over the photon polarizations and averaging upon the J/ip spin, integrating over the entire 
three-photon phase space and multiplying a symmetry factor ^, we are able to determine 
each of the short- distance coefficients tabulated in ([T]). We readily confirm the existing 



0(a s ) and 0(v 2 ) coefficients [lOj, ]15 



We use the built-in function of the mult i- dimensional global adaptive integration algo- 
rithm in Mathematica to compute the one-loop coefficients. We finally nail down the 
value of the desired constant: 



G = 68.913, (3) 

with the fractional error estimated to be smaller than 10~ 4 . 

Since the logarithmic term in the 0(a s v 2 ) coefficient is overwhelmed by the G, from 
now on we will set fif = m c in ([I]). The full NRQCD prediction now is proportional to 
1 — 4.02 + [— 5.32 + 21.9A a s (fj,)]{v 2 ) j/ ^p. It is clear to see that the new piece of correction 
yields a surprisingly significant and positive contribution. With a reasonable choice of fi, it 
has strong potential to counterbalance the known large negative corrections. 

In the numerical analysis, we take the values of the NRQCD matrix elements from 
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Ref. [30j: 

| (0| X V • e"V| J If) | 2 = 0.446 GeV 3 , (v 2 ) m = 0.223, (4) 

which are fitted through the decay J/ip — > e + e~ accurate through the relative order v 2 . We 
take the charm quark pole mass as m c = 1.4 GeV. The Gremm-Kapustin relation then 
implies that the value of (v 2 ) j/$ is consistent with that given in PJ. 

We take the fine structure constant a to be 1/137. There exists some ambiguity in 
choosing the strong coupling constant a s (/i). Presumably, one may think that any scale 
between 2m c /3 and 2m c GeV could be equally acceptable. 

Taking /x = m c in consequently a s {m c ) = 0.388, which is calculated through the two- 
loop renormalization group equation with A^ D 3) = 390 MeV, we then obtain the NRQCD 
predictions for the decay rate of J /if) — > 37 at various level of accuracy. At LO order, the 
prediction is 6.01 eV. After including 0(a s ) correction, the prediction drops to —3.36 eV. 
If further including the 0(v 2 ) correction, we then get —10.48 eV. However, once the new 
0(a s v 2 ) correction is added, we end up with the reasonable value of 0.91 eV. Dividing these 

we 



predicted partial widths by the total width of J/ip, whose latest value is 92.9 keV [31 
then find the Br[J/V -»■ 37] to be 6.46 x 10~ 5 , -3.61 x 10" 5 , -11.28 x 10" 5 , and 0.98 x 10~ 5 , 
respectively. It is amazing that after incorporating the new 0(ouv 2 ) correction, the full 
NRQCD prediction agrees with both the CLEO-c @ and BESIII data quite well. 

One might be curious about the sensitivity of our predictions to the renormalization 
scale /1. To address this question, in Fig. [T] we explicitly illustrate the scale dependence 
of various NRQCD predictions for Br[J/ijj — > 37] in the range between 1 and 3 GeV. For 
reader's convenience, the BESIII measurement 0] is also juxtaposed in the plot. One readily 
observes that the full NRQCD prediction exhibits much flatter /i-dependence than that only 
incorporating the 0(a s ) correction. Finally, it is interesting to remark that, the agreement 
between the full NRQCD prediction and the data appears to be rather stable in a window 
of between 1.1 < // < 1.5 GeV. 

In summary, in this work we investigate the 0(a s v 2 ) correction to the rare decay process 
J /if) — > 37 in the NRQCD framework. Thus far, the large negative first-order radiative 
and relativistic corrections have prevented one from making any sensible prediction for the 
corresponding decay branching fraction. It is very encouraging that this new ingredient of 
correction turns out to be significant and positive. Including this 0(a s v 2 ) correction appears 
to be crucial to reconcile the NRQCD prediction and the recent CLEO-c and BESIII experi- 
ments. The satisfactory agreement between the data and theory may lend some support to 
the values of the NRQCD matrix elements obtained in Ref. [30j. Finally, we remark that 
further studying the 0(v 4 ) and 0(a 2 ) corrections for this process, may turn out to be quite 
enlightening. 
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